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Abstract 

We study the equilibrium Gibbs states for a Boson gas model, defined by Bru and Za- 
grebnov p], which has two phase transitions of the Bose condensation type. The two phase 
transitions correspond to two distinct mechanisms by which these condensations can occur. 
The first {non- conventional) Bose condensation is mediated by a zero-mode interaction term in 
the Hamiltonian. The second is a transition due to saturation quite similar to the conventional 
Bose-Einstein (BE) condensation in the ideal Bose gas. Due to repulsive interaction in non-zero 
modes the model manifests a generalized type 111 i.e. non-extensive BE condensation. Our 
main result is that, as in the ideal Bose gas, the conventional condensation is accompanied by 
a loss of strong equivalence of the canonical and grand canonical ensembles whereas the non- 
conventional one, due to the interaction, does not break the equivalence of ensembles. It is also 
interesting to note that the type of (generalized) condensate, 1, 11, or 111 (in the terminology 



of van den Berg, Lewis and Pule [jTT], |T^, |T5), has no effect on the equivalence of ensembles. 
These results are proved by computing the generating functional of the cyclic representation of 
the Canonical Commutation Relation (CCR) for the corresponding equilibrium Gibbs states. 
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1 Introduction and setup of the problem 



In recent years, the phenomenon of Bose condensation, described first by Einstein in 1925 
has become an active area of research, both experimentally and theoretically. An example is 
the existence of a new kind of condensation which was recently theoretically discovered by an 
analysis of the thermodynamic behaviour of the Bogoliubov Weakly Imperfect Bose Gas 0-0 
or of some specific Bose systems with diagonal interactions P, ^ . This new Bose condensation, 
denoted as non- conventional Bose condensation, is in fact induced by a mechanism of inter- 
action whereas the conventional one, i.e. the Bose-Einstein (BE) condensation, appears by a 
phenomenon of saturation, i.e. by the existence only of a bounded critical density |]T0[-|18|. In 



fact, the Bose condensation occurring in the Huang- Yang-Luttinger model and in the, so-called. 
Full Diagonal Model, studied in great detail in |TP[-|^, should also be considered as examples 



of the non- conventional type, since in the both cases it is due to the interaction in those models. 

The analysis of the effect of the conventional BE condensation on the equilibrium states was 
initially worked out by Araki and Woods in the case of the Perfect Bose Gas (PBG) [Q, and 



further refined in [^-|^. A well-known model that exhibits non- con?;entzona/ condensation is 



the Bogoliubov model 0-0) [fl- As a complete and rigorous analysis of the Gibbs states of the 
Bogoliubov model is beyond the reach of current techniques, we propose to analyze the effect 
of the non- conventional Bose condensation on the Gibbs states in the simpler model defined in 
see ( |1.1|) , in which the both kinds of Bose condensation occur. 

The model we consider is a system of spinless bosons of mass m enclosed in a cubic box 
A C M'^ of volume V = |A| = L'^ centered at the origin with a Hamiltonian of the form 

Hi = Tj, + Ul + Ui = Hl + Ul (1.1) 

with 

Ta = ^ Skaldk, £k = h'^k'^/'2m, for all k 0, 

fceA*\{o} 

Ul = eoa*Qao + ^aSaS^oOo, £o < 0, 5-0 > 0, (1.2) 
= 2^ 9kalalakak, g+> gk> 9-> 0. 

fc6A*\{0} 

The sums run over the set 

A* = IkeR"^: k^ = — n„ = 0, ±1, ±2, a = 1,2, ...,d 

i.e., we consider periodic boundary conditions on dA. We denote the corresponding one-particle 
Hilbert space by (A). Here af = {ak or a^} are the usual boson creation and annihilation 
operators for the one-particle state i/jj. (x) = V~2e^''^, k G A*, x G A, acting on the boson Fock 
space J^f = (L2 (A)) over (A): 

J-f nP, (1.3) 

n=0 
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where 

Hi") = (L^(A"))_ (1.4) 



' symm 

is the symmetrized n-particle Hilbert spaces appropriate for bosons, and TC^^^ = C We denote 

by 

the restriction of an operator A acting on the boson Fock space JF^ to 7i^\ 

This Bose model was introduced and studied in [8]. Clearly, it can be considered as a 
perturbation of the kinetic-energy Ta with diagonal interactions in modes k = {U^) and 

M0(f/1). 

The main interest of this model is that it exhibits two phase transitions accompanied by 
the formation of non- conventional and conventional Bose condensation. The first is due to 
the negative effective excitation energy Eq < 0, which leads to a macroscopic occupation of 
the zero-mode in some interval of negative chemical potentials. The second is a conventional 
condensation due to saturation. Notice that the second repulsive term {go > 0) in prevents 
the Bose gas from collapse, i.e., it keeps the particle density finite. 

We now summarize the main results of IBI , where it was shown in detail that the model Hj^ 



Op displays a two-stage Bose condensation. Let n and 9 = (3 ^ denote the chemical potential 



and temperature, respectively. Furthermore 

fcgA* fceA* 

is the particle-number operator and {—) h{ (Z^' represents the grand-canonical Gibbs state for 
the Hamiltonian Hj^. Define 

(/?, ^) = ^ / - 1) ' dk (1.5) 

pf(/5)^ supp^(/3,/i) = -^ / (e^^^-l)"'rfA;<+oo, rf>2. (1.6) 



and 



and pf are the density and the critical density of the PBG, respectively. The following 
results are proved in [8]: 

• The model has well-behaved thermodynamics, i.e., the pressure exists, for the tem- 
perature 9 = I3~^ > 0, and chemical potential /i < 0. We denote this domain by 
Q = {(^, /i) : 6* > 0, < 0}. Notice that the same is valid for any finite eo £ 1^^- 

• There is no condensation for p < but condensation occurs for Sq < p < Q. More 
precisely, one has a macroscopic occupation of the A; = mode, given by 

pUP^p)^ lim^^^^^(A/i)=max|o,^^|, (1.7) 

i.e., there is Bose condensation due to the instability implied by the negative excitation 
energy < 0, thought of as being induced by an interaction mechanism {non- conventional 
condensation) . 
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For d > 2, the non- conventional Bose condensate density Pq(/?, /i) ( |1.7|) and the total 
particle density 

p'{(3,^^)^ lim^^^^^(/3,/x)=p^(/3,/x)+Po'(A/i) (1-8) 

attain their maxima at /i = 0. For densities exceeding a critical value, 

p>pi iP) ^ sup p) = hm p' iP, p) = pf - £^ < +00, (1.9) 

the Hj^ model (1.1) manifests a generalized type III (i.e., non- extensive) conventional BE 
condensation: 

(ftp) - ito urn 1 5: (;v.>„, ^ { ° ; s |g . (i.io) 

{feGA*: 0<||A:||<<5} ^ 

The term non-extensive refers to the fact that no single mode has a macroscopic oc- 
cupation of particles. For Eq < Q this conventional condensate coexists at p = with 
the non-conventional condensate Pq (/?, p = 0) in the mode A; = 0. Notice that the con- 
ventional BE condensation (1.10) appears in spite of the repulsive interaction (1-2) 
between bosons in modes A; 7^ 0. But it is because of this repulsion that the condensation 
is non-extensive. 

Remark that formula ( p..7|) is also valid for > 0. In this case one gets pg (/3, p) = 0, i.e., 
the model ( |1.1|) manifests, for > 2, only the non-extensive conventional condensation 

(H). 



More details about the non-extensive BE condensation one finds in Appendix A of the present 
paper. These results are an extension of those of |p. They are indispensable for calculation of 
the generating functional for the model ( p..l|) . 

We conclude this introduction with a few remarks. The first concerns the effect of the 
repulsive term f/^ in (1.2). It is known that for (7_ > 0, this interaction converts the conventional 
condensation from type I (macroscopic occupation of bounded number of modes 7^ 0), e.g., 
a single mode = such as occurs in the PBG , eo = fl'o = 5'+ = 0), into one of type III (no 
macroscopic occupation of a single mode, accumulation of a finite fraction of the particles in 
an infinitesimal interval near = 0) 0. The simplest example corresponds to the PBG 



(f^o = fi'o = ^7+ = )in an isotropically dilated container, when the macroscopic occupation of 
the single mode = is transformed by the pure repulsive interaction, go > 0,g^ > 0,Eo = 0, 
into a non-extensive BE condensation [^]. We stress here the isotropic shape of the container, 
since the conventional condensation is so subtle matter that the PBG itself manifests the non- 
extensive BE condensation, if for example, dilated rectangular box has a highly inisotropic 
growth rates for the edges H, [l^, |l^, |T5| . 



Our second remark concerns the dimension dependence of the phase transitions. In con- 
trast to the conventional condensation caused by saturation for d > 2, the non- conventional 
condensation ( |1.7| ) is due to interaction, and it exists for all dimensions, including d = 1, and 2. 
This is another indication that the simplified model (|1.1| ) with diagonal interaction is similar 
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to the Bogoliubov Gas 0-0. Moreover, in contrast to conventional BE condensations, the 
non-conventional condensation may emerge as a first-order phase transitions: the Rose con- 
densate density appears discontinuously, see for example the thermodynamic behaviour of the 
Bogoliubov gas [0; or the Huang- Yang-Luttinger model, see [|I9|, ^ and |2^. This 

transition in the simplified model studied here is continuous however, see (|1.7|) . 

The bulk of the present paper is devoted to the study of the Gibbs states of the model 
In particular we shall calculate the generating functional of the cyclic representation 
of the Canonical Commutation Relations (CCR) for the Gibbs states of the model ( p.. I]) , a 
method introduced in 1963 by Araki and Woods ||2^. From the generating functional it is 
then straightforward the read off properties such as the breakdown of the strong equivalence of 
ensembles as was done in [^]-|^, We shall see that the two phase transitions have their 



distinct effects on the generating functional. 

Before we embark on the actual calculation we present, in Section |^, the relevant known 
properties of the generating functional of the cyclic representation of the CCR for the Gibbs 
states of the PBG [^. In Section |^ we calculate in the thermodynamic limit the grand- 
canonical generating functional for the Gibbs state (— (/?, /i) associated with the model 

(|1.1| ) for a fixed chemical potential < 0, or a fixed density p < pi In the next Section 4 
we determine the generating functional for a fixed particle density p > pi{P) (1-9). In Section 
^ we summarize our conclusions and formulate some tentative generalizations. Some technical 
results are collected in the Appendix A. 



2 Generating functionals 

The purpose of this section is to review the characterization of (Gibbs) states of a Bose system 
by their generating functional, a method originally introduced by Araki and Woods in the case 
of the PBG |^2| . For each Gibbs state there is a representation of the Canonical Commutation 
Relations (CCR) given by the GNS construction. For a complete description see [^, and also 



HI" mi for ^ detailed analysis of the PBG Gibbs state. Here, we only present a quick overview. 

Let M be a complex pre-Hilbert space with the corresponding scalar product (., .)jv/- We 
consider a representation of the CCR over M given by a map h W (h) from M to a space 
U (7i) of unitary operators on a Hilbert space Ti satisfying 

W {hi) W {h2) = exp Im {h, /i2)m| W {h + , (2.1) 

and such that the map \ ^ W {\h) from M to t/ (H) is strongly continuous. By Stone's theorem 
[27], the continuity implies the existence of self-adjoint operators R{h) such that 

W{h) =exp{iR{h)}. (2.2) 

The R{h) are called the field operators and can be interpreted as the random variables of a 
non- commutative probability theory, since by (2.1) one gets 

[R{hi),Rih2)]=ilmih,h2)^,. (2.3) 
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Note that the map h ^ R (h) is a hnear over R. For h G M, we can now define the creation 
and annihilation operators a* {h) and a {h) = (a* {h))* by 



a* (h) 



1 



{R{h) -iR{ih)} , a{h) 



1 



{R{h)+iR {ih)}. 



(2.4) 



V2 



V2 



A representation of the CCR is called cyclic if there is a vector Q in H such that the set 
{W (h) is dense in Ti.. Such Q is called a cyclic vector. It can be shown that, for every 

regular Gibbs state (■), there is unique (up to unitary equivalence) representation of the CCR 
with cyclic vector Q such that 



The generating functional plays the same role for a state on the CCR algebra, as the charac- 
teristic function for probability distribution, see [27]. 

Theorem 2.1 (Araki-Segal) Let E he the generating functional of a cyclic representation of 
the CCR over M. Then it satisfies : 



{Hi) for h e M, the map A — > E {\h) from C to R is continuous. 

Conversely, any generating functional E : M ^ C satisfying {i), (ii) and (Hi) is a generating 
functional of a cyclic representation of the CCR. 

Our concrete setup will be as follows. For a (sufficiently regular) finite volume, A C R*^, the 
grand canonical Gibbs state (■)^(/9, /i), is defined on the set of bounded operators acting on 
the boson Fock space JF^ = JF^ (L^ (A)) over (A), see (1.3). In order to analyze the state 
(■)^ {P, ii), we use the Fock representation W^-^a of the CCR over the pre-Hilbert space M = V\ 
(the space of the C^(A)-functions with compact supports contained in A). Its generating 
functional (2.5) is equal to E^s(/i) = e~4ll'*ll , where cyclic vector Q is vacuum in TC — J-^'- 

a{h)Q, — for any h e V\. Since I>a is dense in (A), one can extend W-^^ to the later. We 
shall calculate the generating functional 



(exp{ii?(/i)}) = (Q,iy(/i)Q)^. 



The generating functional of the representation is defined by 



E {h) = {Q, W (h) n)^, he M. 



(2.5) 




l,s=l 





(2.6) 



and study its thermodynamic hmit (A | R'^) . 
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3 Gibbs state and non-conventional condensation 

Recall from Section that condensation in the exactly solvable model Hj^ ( |1.1D occurs in two 
stages: for intermediate densities p < pi (P), i.e., for negative chemical potentials Eq < p < 0, 
one has only non-conventional Bose condensation in the k = mode due to the diagonal 
perturbation (1.2) of the PBG (cf. (p..7|)), whereas for large densities p > pi (/3) {p = 0), this 
condensate coexists with conventional (type III) generalized BE condensation corresponding to 
the standard mechanism of saturation, see (1.9) and ( |1.10| ). 

In this section we study the influence on the corresponding Gibbs state of the first stage 
of condensation: the non- conventional one ( |1.7| ) which appears for a fixed chemical potential 
£o < P ^ 0. Following |2^-|2^, we use the Fock representations of the CCR over the space 



V\ of C°°-smooth functions with compact support contained in A (Section ^ and we define by 

EiiP,p;h)^(w^S(h)) ^iP,p),heVA, (3.1) 

the grand-canonical generating functional of the model (1.1): 

Hi=Yl {^'^"^^'^ + 1^^^"'} =Y1 ^k, ^Mo = h''e/2m > 0, £o < 0. (3.2) 



feGA* keA* 



Here the operators W (h) for h G Pa? are defined by (2.2)-l 

Note that the boson Fock space J-"^ (1.3) is isomorphic to the tensor product: JF^ ^ (g) JF^, 

keA* 

where = J^^ {Hk) is the boson Fock space constructed on the one-dimensional Hilbert space 

^^ = {^^^'^}Aec- (3-3) 

Then using the Fourier decomposition 

/I . 1 
dxh {x) a* {x) = ^Yl (^'''"' ^)l2(a) i'^k) = ^Y1 ^^'^k 
^ k&A* ^.g^. 

dxh{x)a{x) = -j=jY^ (/;,,e^'=^)^2 a(^J = — = ^ Mfc, 
we can write the generating functional Ej^ (/5, /i; h) (3.1) in the following form: 



(3.4) 



keA* 



n 



° fceA*\{o} 



JJ (e^y^^'''=^^'=''*^'^ . (3.5) 



Next, we study the two factors corresponding to cases = and A; 7^ separately. Denote 
by P = UacR'* space of C°°-smooth functions on M'^ having compact support, and by 



Theorem 3.1 Let Eq G M} and Qq > 0. Suppose that < Q- < gt < 9+ for k G A*\ {0}. Then 
for fi < and any h in the space V one gets that: 
{i) for the mode k = 



lim (e^ 



hoao+hoa'^ 



he 



(3.6) 



where the non- conventional Bose- condensate density pg (Z^,/^) is defined by (1.7); 
(ii) for the second factor in ( \3.3{ ) we have 



lim Yl ('e^*^''*""^^"'^' 



= exp I \\hf - ^Ap^^ {h, h) 



/cGA*\{0} 

where the sesquilinear form A^^^j (m, v), for u,v eV, is defined by 



(3.7) 



{27r)'^ J e^(^fc-^) - 1 



dk. 



(31 



Proof, (i) Let {V'n}n>o orthonormal base of eigenvectors of the operator agao: 

Then one gets 



X = Tr^B ( e'^i^O-I^No) ^-^{hoao+hoal 



+ 00 



n=0 



By the Baker-Campbeh-Hausdorff formula: 
we obtain that 



(3.9) 
(3.10) 



I 

exp <^ — = (/looo + hoal) 



'2V 

Therefore, since 



h^al \ exp <^ — ==/;,oao \ ■ (3.11) 



/2V 



by (3.11) the trace (3.9) equals: 



n=0 

+cx> 



I /in 



n=0 



«=0 



2V / (/!)'(n-/)! 



(3.12) 
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Using the Laguerre polynomials 



Lniz)= 7^ — '■ (—zy, n>0, 



(|3.12| ) can be rewritten to give 



n=0 



2V 



Consequently, we obtain 



2n V ' ^ 



n=0 

Notice that the probability distributions in (3.13): 



^ , , 0<n/V<x 

^v(x) = -^3^ (3.14) 



^^-/3y[(.o-/.-#)f+f (f) 

n=0 



satisfies the Laplace large deviation principle [28, 



A 



lim / Fv{dx)f{x) = f{max{0,- — —}), (3.15) 

go 



for any bounded continuous function / on M^. 
The Laguerre polynomials have the property 



oo ^ 

Jim^ {z/n) = Jo {2^) = tt^ {-^f , ior z e C, (3.16) 

°° 1=0 ( 



as entire analytic functions in C. Here Jo (x) is the Bessel function of order 0. Using this and 
(3.15), we find the thermodynamic limit of (|3.13| ) to be 



lim ^e^(^°'^»+'^°"S))^^ = Jo I d2piie,f^) 



2 

ho 



for any h eV, since ho = ho for A sufficiently large. Thus, by definition ( pTTl) , we deduce (|3l6| 
{a) Similar to the proof of ( p.l3| ) we get for k G A*\ {0} that 

■K^^^-/^v^f(e.-^-#)f+f (f; 



n=0 
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Let 



7fe = e 



+ r , 



Efe-M-fp')"] 



(3.18) 



Since 



we readily get that 



+00 



z] s 



n=0 



n=0 



l-S 



exp < ~z 



1-s ' 



{1 . 1 

see (3.7). 

To show that one gets the same hmit for hm J| we define 



(3.19) 



(3.20) 



feeA*\{o} 



Tk{tk) = e 



+00 



2V 



-pv[{ek-H-^)^ + '-^{^? 



ge-/3v[(.,-M-#)f+^(fP] 

n=0 



(3.21) 



Therefore, Tk{tk = 0) = 7^ and 1^(4 = g^) = T^. Since Tk{tk) G C°°(M^) for each k e A*\ {0} , 
then to prove that hm Y{ coincides with (3.20) it is sufficient to estimate the 

^ fcGA*\{0} 

asymptotic behaviour of derivative 

dt,Tk{tk = 0) 



Afe(y) 



for y — > 00. Here 

+00 

n=0 
+00 

+ E 

n=0 



{Ln 
{Ln 



\hk 



2V 



|g,-/3[(..-.-#)n]| 
U=0 J 



ne-/3[(^.-/'-#)-] (-^ 

^ \ 2V 



2^ , 



(3.22) 



n=0 



\hk 



2V 



- l}e~^[(^^~^~^)"] 



E 

n=0 



2V^ 



.(3.23) 



Since (3.16) imphes the convergence of derivatives, one gets the estimate 

< C^o, z e [0,2:0], 



for any ^0 > and 11 <m. Therefore, in this domain we have : 



Ln{-) 

m 



< CzoU—. 

m 



(3.24) 



(3.25) 



Let no{V) = [V^~^~\ for some S G (0, 1). Here [x] denotes the integer part of the real x. Then, 
since gk < g+, by virtue of (3.25) we can find Ck > such that for any /x < one gets the 
estimates : 



no(V) 

E 

n=0 



{Ln 



\hi 



2V 



l}e-^i( 



+ 00 r , 



n=0 



<c, 



4^2 



(3.26) 



and 



no{V) 

E 

n=0 



I /it 



2V 



l}e-^[( 



9fe ' 
2V . 



+ 00 

E 

n=0 



2V . 



2V 



2\ n 



< C^.^^[5,/(i/ = Pie, - md'yfiy = Pie, - /i))], 



(3.27) 



where fiPie, - /i)) = (1 - e-^(^'=-^))-i. 

On the other hand, for large n the Laguerre polynomials have the following asymptotics : 



Lnix) 



(n'^nx] 



cos [2Vra - 7r/4] + 0(n"^/^), 



for X > 



Lnix) = e"/Vo(2v/(n + l/2)x + 0(n-=^/^), 

for a < X < 6, a > 0, and 

Lnix) = 1 — nx + O ((nx)^) , 
for nx — > 0. Therefore, we get, for some Diik), D2ik) > 0, the estimates : 



(3.28) 

(3.29) 
(3.30) 



+ 00 

E 

n>no{V) 



{Lr. 



\hi 



2V 



l}e-^[( 



)n] ( ^ 

2V 



r / 

J2 e-^[(^^- 



9k 



n=0 



(3.3i; 



and 



E 

n>no(y) 

< /^2(A;) 



{Ln 

JL 



2y 



l}e-^[( 



E 

n=0 



-/3[(efe-M-|^ 



-/3(..-M)no(V)52^(^^^(^^_^)), 



(3.32) 



Since < 0, then for large V the estimates (3.31) and ( p.32| ) are of the order 0(6"^'-^'="^^^*^ *') 
for some < 5 < 1. Taking into account ( p. 261) and (|3.27|) one concludes that AkiV), and 
consequently (F^ — 7^), have for large V the order Oie~^^''V~'^). This implies that for any 
heV 

hm n rfc=lim n 7fc(l + 0(e-'^^'=r-2)) =lim J] (3.33) 



A:6A*\{0} fceA*\{0} 

which, by (3.20), proves the assertion 



fceA*\{o} 
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Remark 3.2 The conditions < g- < gk ^ 9+ for k G A*\ {0} can be relaxed. If < gk = 
gk (V) < -fkV""' for k e A*\ {0} , with < a+ < 1 and < 7^ < 7+, then still holds. 



The proof is obtained by following the same line of reasoning as in the proof of Theorem 3. 1 



Remark 3.3 The first result of the Theorem \3. j| , i.e. ^3. 6i ), is similar to the result for the 
PBG at densities p > Pc (P) (1-6) in the canonical ensemble {P,p) (see j^, [E^/j; 

lim (e^('^-»+^°'^S) \ ^) ^Ai^ ^ = Jo ( J2p^{f3,p) ho ) 

A \ /ta a ^^^^^^ (^{e~^^A}(")j / 

where n = [Vp] is the integer part ofVp and Ti.^^^ is defined byl \1.4\ ). Here 

Pro ^ P.m r /«S«o\ , , lTr^(;o({aSaoe-/^^A}^"^ 
^Ffa„\ „F r a\ / \ f ^ ^hni ^ . 



is i/ie canonical density of conventional BE condensate in the PBG. It should be stressed, 
however, that the limit (^.(\ ) is computed in the grand- canonical ensemble (/?, p) for the non- 



PBG 



Corollary 3.4 Let < gk = gk {V) < 7^1/"* for k e A*\ {0} , with ak < a+ < 1 and 
< 7^, < 7^. Then for Eq E M.^ , p < 0, and h eV, one has 



(/?, p; h) = lim p- h) = Jq (^^2p^ (A/i) |/io|^ exp \-\ ||/if - ^A^,^ (/i, /i) 



4 " " 2 



(3.34) 



where p^ {f3,p) and Ajs^^ {hi, h2) are respectively the non-conventional Bose condensate density 
( \l.Ti ) and the positive closable sesquilinear form (3.8) with domain V. 



Proof. See Theorem |3.1| and Remark p.2| . 

The last result is obtained for any fixed chemical potential p < 0. It can be shown that 
for finite volume and any density p > there is a one-to-one correspondence between p and 
chemical potential pj^ {l3,p), which is solution of the equation 

pi(Ap) = (^)^/Ap)=p. (3.35) 

From Corollary 3.4 we have an explicit calculation of the grand-canonical generating functional 
Ej^ {(3,p; h) in the thermodynamic limit. However, for a fixed total particle density p > one 
has to evaluate the following thermodynamic limit 

E^ (/3, p; h) = lim E^ pi (/3, p);h), heV. (3.36) 
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This is not guaranteed to equal the hmit ( p.34|) with chemical potential fi^ p) =1™ /^i (/^' P)^ 
since the map p ^ fi^ (/?, p) fails to be injective: 

/ {13, p) = hm pi p) = I < U < llJl , (3.37) 

because of conventional BE condensation for dimensions d > 2. In this case the total particle 
density p^ {P,p) (1.8) is saturated for p = 0: there is a finite critical density of particles 
pI{9) , ci. (1.6), (1.9). The question of equality of the hmits ( p.34| ) and ( p.36| ) will be considered 
in the following paragraphs. 

In fact it is linked to another question, which concerns the relation between thermodynamic 
limit of the canonical generating functional defined for the total particle density p > by 

K,can P; h) ^ ^- heVj,,n= [Vp] , (3.38) 



and the grand- canonical generating functional W {(3,p;h) defined by (3.36). In other words, 
in the thermodynamic limit, for the same particle density p, the canonical ensemble may not 
yield the same equilibrium state as the grand- canonical ensemble, i.e. one may have 

(/?, p; h) =\im Ei^,,„ (/?, p; h) ^ (/5, p-h),heV. 
To answer to these questions we notice that 

e-/^^[-Mi(4f )]Ei,_ (Af;/.) 
Ei pi (p) ; h) = — ^ ^ ^ , n = [Vp] , (3.39) 



+ 00 

E 

n=0 

where fl (/3,p) is the free-energy density associated with the Hamiltonian Hj^ (1-1)^ 



e-^""^ \ , p > 0, n = [Vp] . (3.40) 
It is known from [8] that : 

(^) {/a (/^) p)}acm<* the family of strictly convex functions of p > 0, and this is also valid for 

/^(/?,p)^ lim/l(/?,p), (3.41) 

A 

the free-energy density fj^ {(3, p) ( 3.40) in the thermodynamic limit, but in the smaller domain: 

p<pi{P) ; 

(b) the convergence ( |3.41|) implies in this domain (Griffiths lemma, see e.g. [§]) the limit : 

dpf (/3, p) = hm dpfi (/3, p) ^ pL P) < 0, (3.42) 
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which coincides with ( p.37|) : 

/^L(Ap)=/(Ap), (3.43) 
for p < p[ , where this function is one-to-one ; 

(c) the grand-canonical pressure, {P,p) =hmpy^ {P,p), is the Legendre transform: 

p'{P,p)=snp{fip-f{P,p)} = {fip-fiP,p)} , (3.44) 



p>0 

and 



f{P,p)=snp{pp-p'{P,p)} = {pp-p'{P,p)} (3.45) 



M<0 



f^=p'W,p) 



where p^ (P,p < 0) < pi (/3) is the function inverse to the injection (3.43): 

p' (/3,/(Ap)) =p; 
(d) for p> pI (/?), see ( |LTO|) and ( p37D , the hmit { ^AJ\ j is equal to 

f {(3, p) =sup {pp-p' {13, p) ] = V (/3, = 0) , (3.46) 

i.e. the free-energy density is not a strictly convex function in this domain; respectively, 

p^{P,p = 0) =sup {—f^{P,p)\ = —f^{l3,p>pi{P)), which means that the pressure and 

p>o 

the free-energy density are always related by the Legendre transform: weak equivalence of 
ensembles. 

By virtue of (a)-(c) we can now apply the Laplace large deviation principle [28, 29| to 
calculate the limit of ( p.39| ) in domain p < pl{(3): 

iP, p- h) = / (/3, p)-h)= (/5, p; h) . (3.47) 

Notice that for d = 1,2 one has p^ = +cxd, see (1.9). Therefore, Corollary together with 
( |3.47| ) imply the following result. 

Theorem 3.5 Let p < pi (P). Then 

{p,p<piiP);h) = Jo U2piie,p'iP,p)) ho ) exp \\hf - ih,h)\ , 



. 4 " " 2 

(3.48) 

for d > 1, eo e R^, < gk < ik^^ for k e A*\ {0} , with < a+ < 1 and < 7^ < 7+ 
and h E v. Here A^^^ {hi, h2) and p' p < pi (/3)) < are defined respectively by (3.8) and 

Consequently, the equality (|3.47|) shows the strong equivalence between the canonical en- 
semble (/5, p) and the grand-canonical ensemble {[3,p) for p < pi{(5) (i.e. p < 0): in the Hj^ 
model (|1.1|) for a fixed total particle density p < pi{l3) the Gibbs state in the grand-canonical 
ensemble coincides with the one in the canonical ensemble. 
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However, contrary to the non- conventional Bose condensation (|1.71 ) the conventional BE 
condensation 'pq{i3,p) > ( |1.10| ) violates this strong equivalence. Indeed, by virtue of (d), 
see (|3.46| ), the limiting measure (|3.39|) relating two generating functional is not degenerate in 
domain p > (/?). Therefore, similar to the PBG the existence of the critical density 

pI (P) implies that for p > p[ one has 

E^(Ap;/i)^EL„(/3,p;/i). (3.49) 
In the next Section 4 we show that, in contrast to (|3.47| ), for p > p^ {f3) one also gets 



(3.50) 



4 Gibbs states and conventional condensation of type III 



Since p^ (/5) < +oo (1.9) only for c/ > 2, we consider A C M'^^^. In the interest of simplicity 
we restrict ourselves to a cubic box of the volume V = |A| = L'^, and we put gk = 9 > 0, 



G A*. Notice that our reasoning in the proofs of Theorems |3.1| and p.5| used that /i < 0, and 
that //j^ (/?, p < Pc (/?)) < 0, for large V. For p > pi (/5) it is not the case, see [Appendix A . 
This difference modifies essentially the calculations of thermodynamic limit of the generating 
functional. 

1°. Our first step is to refine, for V +oo, the asymptotics of the chemical potential 
A^A i.PyP)y which is solution of equation (|3.35| ): 



pi 



(/3,pi(/?,p)) =:i J2 {Nk)Hi{P,pi{P,p)) + 



No 
V 



(/5,pi(/3,p)) =p. 



Hi 



fceA'\{o} 

For a strictly positive g > this is done in [Appendix A| , see Theorem [A. 7] : 



(4.1) 



where the constant C = {2m/fpf^^ / [g 2'^-^7T'^/^d{d + 2)T{d/2)] > 0, and T{z) is the Euler 
gamma function. Hence, the chemical potential /^j^ (/3,p > pi (/?)), is non-negative for large V. 
By virtue of (3.13) and ( p.l7[) this observation motivates to represent the generating functional 
([3.5[) in the form : 



K {P,pi {P,p);h) 



/ Hn 



X 



n 



g V2V ^ 



l(A) 



^''^'^-^^^^l)) (/?,pi(/?,p)) 



X n (e^^'^"'''^'^^^)^,(A/^i(Ap)), 



(4.2) 
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with 



^ |fc e A*\{0} : su-l^i{P,p) - t < 0} 



^[ke A*\{0} : e,-^^i[p^p) - f > 0} 



(4.3) 



Remark 4.1 As it is shown in \Appendix Al for p > {j3) and Qk ^ 9 > the non- extensive 
condensation pg {f3,p) is concentrated on the set d[!^K- 



~pi{P,p)= limi Yl {Nk)Hi{/3,pii/3,p))=p-pi{^)>0, 



(4.4) 



(A) 



see Lemma A. 6 in {Appendix A 



2°. Since in the proof of Theorem |3.1| (i) the sign of p is irrelevant, the same hne of reasoning 
gives that for p = pjt^(^P, p > p^ (/?)) , ( |4.1| ), and for any h E V , Eq & M}: 



hm (e^ 



hoao+hoag 



^{(3,pi{f3,p)) = Jo(j2pi{(3,0) 



ho 



(4.5) 



Here the non-conventional Bose-condensate density pg {9, 0) is defined for p = by (1.7). 
3°. Now, using the asymptotics (^4.1| ) we can evaluate the thermodynamic limit of 



n 



e 



(A) 



(4.6) 



see (4.2). Indeed, by inspection of the line of reasoning ( |3.17| )- p.32| ) we find that it is only the 
inequality Sk — pj^ (/5, p) — |^ > which one needs that the limit ( p.33| ) be valid. Therefore, 
taking into account (|4.1|) we get 



hm H 



huak+hkO.!] 



(A) 



H 



^(/3,pi(/3,p)) 



exp <i \\h\f - ^Ap^o {h,h) 



(4.7) 



for /i e D and p > pi (9), with the sesquilinear form A/3,^=o (^i, h2) defined by (3.8). 

4°. Finally, using the asymptotics ( ^.1| ) we have to compute the thermodynamic limit of the 



last factor in ( 4.2 



hm n 



k&D 



(A) 



(A/ii (Ap)). 



(4i 



Because of non-negative p^ (/5,p), and of the non-extensive BE condensation pg (6*, p) (4.4), 
that spreads over the modes k G ( Remark 4.1 ), this calculation is a more subtle matter 
than the thermodynamic limit of ([4.6|) . In particular, the exact knowledge of the asymptotics 
(|4.1| ) of p]^ (/3, p) becomes essential to find the limit of (|4.8|) . 

By (|3T3) we get: 



+ 00 



fhkak+hkal^ 



\hk 



n=0 



2V 



(4.9) 
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where 



UA,k iP, P ; n) 



g-/3[(efc-/.i(/3,p)-#)n+|^n2] 



(4.10) 



n=0 



is the family of probabihty measures {I'A^k {P: P]^)}AcR'KkeA* > cf. (3.14) and (A. 17). Conse- 
quently 



In Y[ rk=(-^j Yl In ^ '^A.fc (/?, P ; ^) ( ) • 

Since /xj^ p > piiP)) 0, the thermodynamic limit of the first term in (4.11) is 

1\ 1 



(4.11) 



lim I I =lim 



4 y (27r) 



(4.12) 



{fc:£fe<5} 

for any h G V. By virtue of the asymptotics (4.1) and by definition of domain D^^^ we get for 
the limit of the second term : 

hm In ^ z^A,fc {l3, p\n)Lr^ 



k&D 



(A) n=0 



2V 



h 



k(s)/Vt/2 



2V 



limj^ In ^ivMs) (V'-'^/3,p;n/r-^) L„ 

se5fl n=0 

Here 7 = 2/(rf + 2), the set : 

Sb^[s = e {0} : £-(2vr)^ {s./V'"'"'f < B 



(4.13) 



(4.14) 



a=l 



and we put k{s) = 2713 /V^^'^' "'1"^. Since for d > 2 one has 1 — 27 > 0, the family of the scaled 
probability measures: 

{vAMs) (^'~'"/5,P;n/V'"")}AcR^.e5B, (4.15) 
see (4.10), verifies the Laplace large deviation principle |2^ with the support at the point 
lim n{y)/V^~"' = {B — ek{s))/g , cf. Theorem [A.6| . This remark, together with the asymptotics 

(3.30) of the Laguerre polynomials and the continuity of h &T>, gives : 



h 



k{s)/V''/^ I 



lim^ InL^ 



h 



{V)=V'-'{B-ek{s))/g 



fc{s)/yT/2| 



sGSb 



-limV V 

1 - 2 



— ek{s)) |^fc(s)/v^/2 
21^" 



2\/ 

|2 



ho 



2 1 



(2vr) 



{k:ek<B} 



{p-p[m- 



(4.16) 
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Here we used that 7 = 1 — 6/7/2 to obtain in the hmit the integral, and ( A.25| ) to get the last 
equality. Taking ( [4.12| ) and ( |4.16| ) into account we finally get for ( |4.8| ): 



hm H 



k&D 



exp 



(A) 



(/5,P)) = exp 



1 h ' 



(p - pi m) 



where {(3,p) is density of the type III BE condensation, see ( p. .101) . 
The results of 1° — 4° can be summed up as follows: 

Theorem 4.2 Let p > pi (P). Then 

E'iP,p;h) ^ \imEi{p,pi{P, p);h) 

ho 



Jo{^J2piiP,0) 
1 



exp 



1„.,,2 1 



2^/3,0 {h, h) 



X exp 



(P - Pc 



(4.17) 



(4.18) 



for d > 2, Eo e M^, gk = g > and h E V. Here Ai^^^ {hi, h2) and pj^ (^(3, p > p^ (/?)) > are 
defined respectively by (3.8) and 



Corollary 4.3 Comparing ( \3. 34 ) and ( U-i^ ) one gets that the presence of the conventional BE 
condensation pg (/5, p) > 0, see ( \1 ■!(]) , ( [i-4 )> implies : 



iP, p; /i) = (/?, 0; /i) exp I - 1 I /Jo I ' (/?, p) | 



(4.19) 



cf. (3.50). This means that at the point p =lim pj^ (/?, p > pi (/?)) = the grand- canonical 

equilibrium state is not unique. There is a family of states enumerated by the BE condensate 
density pg (/5, p) = p - pi (p) . 



This effect is well-known in the PBG ^-^\ 



Proposition 4.4 For isotropic dilation of a rectangular container A C M"'^^ the grand-canonical 
generating functional 



iP,p; h) =lim 



Tr^s ( n e^*^''*'''+'''''^)e"^(^^"''A(/5:P)^A)A 
^ VfceA* / 



Tr^fl ( e-K^A-M^{/3,P)^A; 



lim exp <j -i II /if - ^A^,^p(/3,p) h) j> lim (eV2v 



(/5,PnAp)) 



exp <; -] \\hf - ^A^^^P(^^p) {h, h) \ exp | ~ 



ho 



(4.20) 
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Here fi^ {f3, p) < is solution of the equation 



V 



1 Tr^s (iVAe- 



■l3{TA-fiNA) 



) 



Ta 



Similar to ( \3. 3'1\ ) and ( \4-4\ ), one has: 



and t/ie i?/j»e / conventional BE condensation in the single mode k = 0: 



Po iP^ P) =lip 



\ ^ Ita 



!orp<pl{(5) 
p-p^AP) forp>pAP) 



Hence, similar to (4.19), one has 



(4.21) 



E"' p > pf (/5) ; /i) = (/5, 0; /i) exp -- ^ p^ (/5, p) 



Moreover, the remarks (a)-(d) are vahd for PEG. For p < pf (jS) there is a strong equivalence 
of ensembles expressed by : 



if], p; h) = E^ (/5, p^ {(3, p);h)= E^ (/?, P; h) 



(4.22) 



cf.(3.47). Whereas for p > pf (/3), the functional E^„(/3, p;/i) and the non- degenerate mea- 
sure in ( p.39| ), relating the canonical and the grand-canonical generating functionals, can be 
calculated explicitly p3[-p^. This gives : 



E^ if3, p-h)= irp,pp(^)(t^x) Kan a;; h) . 



(4.23) 



Here 



EL X > pf iP) ; /,) = exp <j -i II /if - 1a^,o {h, h) \ J, (^^2{x-p^m 



and 



for X < pf (/?) 

{p - (/3)}-^ exp ^ - ~ f ^ I do; forx>pf(/3) ' 

P-Pc(/5))J 



(4.24) 



(4.25) 



is known as the Kac measure for p > p^ (/3). Since by (a)-(c) the limiting measure in (3.39) is 
degenerate : 

Kp,pP{/3){^x) = 6{x- p) dx, (4.26) 
for p < p^ (P), cf.(4.22), the representation ( [4.23| ) is valid for any p > 0. 
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5 Concluding remarks 

5.1 To answer the question about the role of the type III BE condensation in determining 
the generating functional ( |4.19| ), let us consider instead of (|1 . 1| ) a truncated model with the 
Hamiltonian H^. 

Since Theorem p.l| is valid for = (or = 0), the generating functional for the model 
has for /i < the same form: 

E°(/?,/i;/i) = hmJJ (^e^l'^'^'^+'^'^'^fe)^ ^(/?,/i) 
fceA* 

= Jo (^^2p[j (A/i) \ho\^ exp \\h\f - ^Ap,, {h, , (5.1) 

as for H{, cf. (|3.34] ). Here pfj (/5,/i) = Pq {l3,fi), see ([1.7|). Again, in this domain we have the 
strong equivalence of ensembles ( p. 47] ): 

E° (/?, p; /i) = E° p° (/3, p) ; /i) = E^ (/3, p; /i) , (5.2) 

where p° (^, p) = p^ (/5, p) and pO (/3) = p^ (/3), see (3.37). 

On the other hand, for p > {(3) the model manifests (instead of the type III) the BE 
condensation of the type I, see P). More precisely, for dilatation of a cubic A there is the BE 
condensation in 2d modes : 

IC2d = {keA*: {fcfjti = (0, 0, 0, K = ±27r/FV^, 0, 0), a = 1, 2, 4, 

such that 

hm i (N.^kJho i/3,pl{/3,p> pI m) = ^(P - Pc (5.3) 

and 

lim i (Ar,eA.\{x:,,u{o}}>^o {P,pl{P,p> P° (P))) = 0- (5.4) 
It corresponds to the asymptotics : 



2d 1 

,5(ftp>P»(^))=e..-^,^^j^^-^ + 0(-) (5.5) 



of the solution of the equation 



77 E {Nk)Ho{f3,plif3,p)). (5.6) 



P V 

keA 



Since — p\{l3, p) > for 7^ 0, following the same line of reasoning as in calculation of (4.2) 
we get for : 



E^(/3,p^(/3,p);/^) = (e^(''»--^'^»'^Sn ^(/?,p°(/5,p)) 

n (ew(^'^^-+'^'='^^)\ (/3,p°(/3,p)) 



^ T T / ' (hi.au 

n (e^(^'^'=+'^'=";^))^ {P.M.P)) (5.7) 
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the limit : 



E° p; h) 



limE" (/3,p);/i) 



X exp 



2pg 0) 



exp < — 



2 1 



Apfi {h, h) 









ho 





{p - pI 



which coincides with the generating functional ( [4. 181 ), or ( f4.19| ). 

Theorem 5.1 The generating junctionals for the models and Hj^ are identical : 

¥P {P,p- h)=E' {P,p- h), 

for any P and p. 



(5. 



(5.9) 



5.2 Notice that the relation between grand-canonical and canonical generating functionals 
(4.23) for the PBG can be seen by means of the identity : 



exp(-lAz^) = I je-'/'Jo (v^ 



(5.10) 



for A > 0. It yields the representation of the grand-canonical generating functionals (4.19) (or 
r^)) via the Kac measure: 



(/?, p;h)= K^,,i^(,){dx) {(3, p = 0;h)Joi V'^ix - pi (/?)) 



he 



(5.11) 



Following (4.23) and (4.24) this gives a temptation to identify 



(/?, p = 0;h) Jo (V^ip-pi m 



hr 



Kan (AP; h) 



(5.12) 



for p > p[ (P), with the generating functional E^^^ (/3,p; h) of the model H{ (or H^) in the 
canonical ensemble. But actually we do not know neither E^^^ (/?, p; h), nor E°^^ {P,P', h), for 
p > pi (/9), to check this hypothesis. 

For p < pi {f3) the representation of the grand-canonical generating functional via the 
degenerate Kac measure (4.26) and the canonical generating functional follows directly from 



(1113), or (13). 

5.3 In the present paper we are not concerned with the explicit form of the OCR representations 
corresponding to the generating functional E^ {j3, p; h) for different densities p. Therefore, we 
limit ourselves only by few remarks : 
- In domain p < p{j3, eo) one has: 



E^ {(3,p<pi {(3);h)=exp 



^ - 2^Am^(/3,p)<o (h, h) \ , 



(5.13) 
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where quadratic form i \\h\f + \Ap^^i(^p^p^ ih,h) is closable. It is known that in this case the 
representation is a factor corresponding to the class of quasi-free states, see I^O] and for 
details. 

- By virtue of ( p.lO| ) the W {P, p; h) in domain Eq) < p < p^ (/?) is the (inverse) Laplace 



ho 



transform of a generating functional with the non-closable quadratic form |A 
|74^_^/(^_p) {h, h). Hence, the representation is not a factor. 

- Because of the term — | Pq p), see ([4.19|) , the same remark is valid for representation 
in domain p> p^. {(3). For details of construction of the representation in these cases see ||24 



Appendix A 

The aim of this appendix is to investigate the type III BE condensate in the model ( |1.1| ), and 
to add some new results to what is known since p, The essential point is to obtain the 
asymptotics of p\ (p), which is the solution the of equation 

.i(/5,pi(Ap)) = :^ WHi(A/^i(Ap)) + (v) (A/^i(Ap)), (A.i) 



fceA*\{0} ^ ' 



see (3.35). We recall that by ( [L.7D one has 

lim/^\ (/3,pi(/3,p)) =p^(0,p^(/3,p)), 

where we put p^(/3, p) =lim p{i(3,p). To simplify the arguments, we consider below a cubic 

A 

box A C R'^^^ of the volume \^ = |A| = L'', and (at the very last moment) we put all constants 

{S'fclfeeA* to be equal to > . 

A.I Let Z}^^ and D^_^^ be two sets defined by 



^(A) ^ 1^ e A*\ {0} -.eu-pi (/?, P) - f < 0} 
^ {fc G A*\ {0} -.Sk-Pi (/?, P) - f > 0} 



(A.2) 



cf. (4.3). Then ( |A.1| ) transforms into 



h E (^^)Hi(/5,pi(Ap)) + ^ E WHi(/?,pi(Ap)) + (^) (/5,pi(/5,p)). 



(A.3) 

Notice that, since Ek^o = 0{V~^/'^) for d > 2, the set D^J^^ = 0, if p{ (/5,p) < for large V. 

Lemma A.I /@/ Lei g+ > gk > 9- > for k e A*\ {0}. T/ien /or any Eq G anc? k G 
one /ifls the estimate: 

/or sufficiently large. 
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Proof. By the correlation inequalities for the Gibbs state uj^ (— ) = {—)h{ (^^^ [S EH)- 

M (X* [Hi ,X])>u;i {X*X) In ^TyS' (^-5) 

where X is an observable from the domain of the commutator (/i) , .] with (jj,) = 
H{ — fiN^, we can deduce 

M (a* [Hi ,a,])> (AT,) In f^^l^'] , (A.6) 

for X = a^; (/c 7^ 0). Since for /c 7^ 0, 



S'fc * 2 

from (A.6) one finds for fi = fij^ (p) 



(iVfe) + l^ 



which gives: 

Now the rest of the proof is essentially due to solution for x > of the inequality (A. 7): 

bk<ln^, (A.8) 

X 

where 

x = ui (Nk) ,h=/3(^ek- pi (p) - ^) . (A.9) 
Since for k G D*^^ we have 6^ > 0, the inequality (A.8) implies ( [A.4| ) for sufficiently large A. 



Remark A. 2 /Q/ Because of repulsive interaction Qk ^ 9- > ^ for k G A*, see (1-2), the 
finite-volume pressure pj^ (/3, p) 0/ i/ie model ( fOl j zs a convex function of fj, E such that 
limAPA (/^'/^ > 0) = +CXO. This means that (in contrast to the PBG) for large p the solution 
A*A (p) equation \A.^ ) should he positive. By convexity : 

< <9^Pa (P, /^) = Pa (P, /^) (A. 10) 

A* 

for p > 0, one (^efo i/iai lim^pj^ {P,^ > 0) = +00. Therefore, \im\fij^ (p) < 0, i.e., the positive 
solution of (A.l) must go to zero in the thermodynamic limit. 

Corollary A. 3 For p > p[ (/3) and V sufficiently large the solution p^ (/?, p > pi {13)) > 0. 
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Proof. Suppose that /x^ (/?,p) < for any p and A. Then dL^^ = 0, i.e., D^^^ = A*\ {0}, and 
by Lemma A.l one gets the estimate: 



fcGA*\{0} 

l_ 1 , /No 



Hi 



V 



Hi 



Then, by virtue of pj^ {l3,p) < 0, we get the following inequalities in the thermodynamic limit: 



limpA f^i P)) + 1™ ^ Y 

Pf (/?) + po' (A p = 0) = 



{P,pi{P,p)) < 



(A.ll) 



see (|1.7|)-(1.9). But this is impossible for p > p{ that proves the assertion. 
Therefore, Remark |A.2| and Corollary |A.3| state that 

hmpi (Ap>p^(/?)) =0. 



(A. 12) 



Corollary A. 4 The representation together with arguments of Corollary \A.^ and (\A.12J , 

allow to refine the localisation of the non-extensive BE condensation 



Po'(/?,p)=liml WHi(Api(Ap)) 



p-^T^ E (iv.)Hi(Api(/3>p))-ir(^) (ApKAp)) 



= P-Pc(/5)- (A. 13) 

A. 2 Our next step is to calculate the asymptotics of p^ {P,p > piiP)) in (|A.12|) . To this end 
suppose that for V +oo it has the form : 

pi{f3,p>pi(f3))=^ + o{V-'), (A.14) 
with 7 > and B > that should be defined from equation (A.l). 

Remark A. 5 Suppose that > 2/d. Since Sk^o = 0{V-^l'^) , then the set D^^^ = 0, for large 
V. Therefore, the same line of reasoning as in Corollary \A.^ produces a contradiction to our 
main assumption : p > pi{P)- Hence, we must have : 

7<2/rf. (A.15) 
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By virtue of additive structure of the Hamiltonian (|3.2|), for any G A* we get that 



AT \ / N \ 1 

{p,^{(p,p)) = /^\ (p,p{^p,p))= J2 nu^,UP,P;n), (A.16) 



with probabihty measures : {i^a.A: P ; ^)} 



AcK'',fceA*- 
g-/3[(.fe-Mi(ftp)-#)n+|^„2] 

i/A,fc (/?, p-n)^ — ^ -. (A.17) 



^ ^-/3[(.,-^^(/3,p)-|^)„+|t„2 



n=0 



From (A.17) it is clear that we have to distinguish two domains : 

^ {A; G A*\ {0} ■.ek-^^i i/3, p) - < O} , 
= {fc G A*\{0} : -/ii(/3,p) - f > 0} , 

cf. (4.3) and (Q. Since L>i^^ C dL'^^ our next statement makes the locahsation of the 
non-extensive BE condensation more precise, cf. ( |A.13D . 

Lemma A. 6 For p > pi (/?) one has : 

Po^ (A p) =hm i Yl {Nk)Hi{/3,f^ii/3,p))=p-piif3)>0. (A.19) 



Proof. By (A. 13) it is sufficient to prove that 

^TV E {Nk)Hi{P.^^iiP,P))=0. (A.20) 



Since < Q- < Qk < g+, hy (A.16) and (A.17) we get that 

(iV.)^i (A pi (A P)) < T;;^^^^^^ < O [V^'') (A.21) 

e-/3Sfcn2/2y ^ g-/3(sfcn/2y+gfcn2/2y) 
n=0 ra=0 

for A; G D^_/^\d^_^'^ and large On the other hand, the set D^_^\d^_^^ has the same number of 
elements as the set 

{^ = K}Liez^{o} : |^ + pi(Ap)<|^^E^' ^y + /^i(^'/^)}- (A-22) 



Since the volume of the elementary cell of the dual lattice A* is equal to {27iY/V, the number 
of points (A. 22) for large V is finite. Together with the estimate (|A.21| ) this gives ( |A.20|) . ■ 
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Theorem A. 7 Let = g > Q. If p > p[{l3), then asymptotics of the chemical potential 
A^A P) ^'^^ ^he form ^A.14\ ) 

7 = 2/{d + 2) and B{(3, p) = ( \ . (A.23) 



Here C = {2m/ffr' / [g 2'^-^7r^/^d{d + 2)r(d/2)] > 0, and T{ z) is the Euler gamma function. 

Proof. One has to tune the values of -B > and 7 > in such a way that to satisfy the equation 
(A. 19) for V —>■ +00. Since we have 7 < 2/d (Remark |A.5|) , by using ( |A.16|) and (A. 17) we get: 

S W^,J(3,pii/3,p)) = 



n=0 



hoo 

n 



I 5Z Yl ^A,fc=2-/vv.-./2 (V'-^^P, p ; , (A.24) 



sG5g n=0 



where Sb ^ {s = {sX=i ^ {0} : £(27r)2 ^ (sjyi/'^-^/^)^ < B} , see (gll). 



a=l 



The sum in (|A.24| ) over 5^ is nothing but the Darboux-Riemann sum converging to the 
integral, when V —>■ +00. Therefore, to get a nontrivial limit in ( |A.24 ) we must choose the 
value of 7 from the condition: 1 — 7 = d'y/2, i.e., ■j = 2/{d + 2) in ( |A.14D . For this value of 7 
and d > 2 one has 1 — 27 > 0. Then the family of the scaled probability measures: 



cf.(A.17), verifies the Laplace large deviation principle [^,|29| . Hence, by a diagonal limit 
involving in ( |A.24| ) the sequence of Darboux-Riemann sums and probability measures, and by 
(|A.19| ), we deduce the equation: 



{d+2)/2 



P-Pi{p) = -^, [ d'k?-^ = { {2m/h^Y" I [g 2'~\'/^d{d + 2)T{d/2)] } B 
(27r) J 9k 

{k:st:<B} 

(A.25) 

which defines the value of B = B{f3, p). This finishes the proof of (A.23). ■ 
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